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SUMMARY 


It  Is  demonstrated  that  from  a  knowledge  of  the  thermal  and 
optical  properties  of  each  layer  of  a  two-layer  system,  together 
with  the  amplitude  and  distribution  of  the  energy  input,  it  is 
possible  to  predict  temperature  rises  at  depth  in  the  second  layer 
of  a  two-layer  system.  By  varying  the  optical  properties  of  the 
first  layer  and  observing  the  temperature  rises  at  depth  in  the 
second  layer,  the  experimental  results  can  be  used  to  verify 
mathematical  expressions  for  optical  and  heat  transfer  processes 
of  the  two  layers. 

Such  a  system  was  devised  to  determine  the  thermal  conductivity 
of  intact  living  skin  by  measuring  the  interface  temperature  between 
a  silicone  rubber  patch  and  the  living  skin  when  irradiated.  It  may 
also  be  used  in  evaluations  of  other  two-layer  systems  where  reflect¬ 
ance,  transmittance  and  heat  transfer  properties  are  known  and  must 
be  accounted  for  in  the  mathematical  model. 
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INTRODUCTION 


Burn  effects  of  exposure  to  non-ionizing  radiation,  whether 
of  low-level  or  thermonuclear  intensity,  depend  upon  the  optical 
as  well  as  the  thermal  properties  of  skin  and  whatever  protective 
cover  it  may  have.  Pigmentation  of  the  unprotected  skin  may  great¬ 
ly  influence  the  magnitude  and  depth  of  penetration  of  incident 
radiation  depending  upon  its  wavelength.  It  is  therefore  impor¬ 
tant  to  be  able  to  measure  the  thermal  conductivity  of  intact 
human  skin  and  to  be  able  to  predict,  from  a  knowledge  of  the  op¬ 
tical  and  thermal  properties  of  both  the  skin  and  its  covering, 
the  effect  of  incident  radiation  of  known  wavelength  and  inten¬ 
sity. 


We  propose  to  adapt  a  two- layer  heat  transfer  system  to  a 
system  for  determining  the  thermal  properties  of  one  layer  of 
material  from  a  knowledge  of  the  thermal  properties  of  the  other 
layer  when  the  temperature  rise  on  irradiation  and  intensity  of 
the  energy  are  known  (1).  In  this  way  a  patch  of  suitable  material 
of  known  properties  may  be  laid  upon  the  living  skin  with  inter¬ 
face  or  surface  temperature  rise  being  measured  during  irradiation, 
making  it  possible  to  calculate  the  thermal  conductivity  of  living 
skin.  The  validation  of  the  proposed  system  entails  the  use  of  a 
series  of  silicone  rubber  patches  of  different  optical  properties 
representing  various  degrees  of  pigmentation.  The  patches  are 
combined  with  a  skin  simulant  to  form  a  two-layer  system,  analogous 
to  a  rubber-covered  living  skin  combination,  in  which  mathematical 


expressions  of  heat  transfer  relationships  can  be  verified  ex¬ 
perimentally.  This  paper  describes  the  heat  transfer  process 
and  may  be  used  in  general  applications  involving  two-layer 
systems  as  well  as  in  our  specific  use. 

PROCEDURE  AND  MATERIALS 

The  two-layer  system  is  irradiated  by  a  graphite  imaging 
furnace  (2)  to  validate  experimentally  the  temperature  rises  at 
depth  in  the  second  layer  which  are  determined  mathematically.  Com¬ 
puter  techniques  are  used  to  calculate  the  optical  properties  and 
energy  distribution  of  the  two-layer  system.  An  equation  for  tran¬ 
sient  heat  flow  through  a  two-layer  wall  as  proposed  by  Griffith  and 
Horton  (3)  is  used  to  determine  temperature  rise  at  depth  in  the  sec¬ 
ond  layer  as  it  varies  with  time. 

Each  of  seven  silicone  rubber  patches  designated  W(white), 

W  +  1 ,  W  +  2,  W  +  3,  W  +  4,  W  +  5  and  Black  is  used  as  the  first 
layer.  Silicone  rubber  was  chosen  because  the  thermal  properties 
are  known,  it  can  be  prepared  in  fairly  thin  sheets,  it  does  not 
wrinkle  or  swell  appreciably  when  irradiated  and  it  forms  a  tight 
contact  with  the  backing  material.  The  thickness  of  the  patches 
is  0.0532  cm  +  0.0008cm.  Each  patch  was  pigmented  differently 
to  vary  reflectance  and  transmi ttance  properties.  The  backing 
material  is  the  skin  simulant  device  fabricated  by  the  Naval  Mater¬ 
ial  Laboratory,  N.  Y.  (4).  Optical  properties  were  measured  with 
the  Beckman  DK-2A  Spectrophotometer.  Reflectance  and  transmi ttance 
were  determined  for  each  silicone  rubber  patch.  Transmi ttance  was 
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measured  for  a  slice  of  the  skin  simulant  0.047  cm  thick, 
approximating  the  depth  of  the  thermocouple  beneath  the  surface, 
and  reflectance  in  the  full  thickness  (1cm).  Using  computer  tech¬ 
niques  the  absorptance  is  calculated  at  each  wavelength  for  each  layer: 
A1  (a)  =  1  -  [R1 (a)  +  T1 (A)]  Eq.  (1) 

A2(a)  =  1  -  [R2( A)  +  T2( A)]  Eq.  (2) 

where 

1  refers  to  top  layer,  2  refers  to  base  layer  and 
A(a)  =  Absorptance  at  each  wavelength  interval,  % 

R(a)  =  Reflectance  at  each  wavelength  interval  =  R(a)SC(a),0/ 
T(a)  =  Transmittance  at  each  wavelength  interval,  % 

A  =  Wavelength  of  radiation,  microns 

SC ( a )  =  Setting  correction  at  each  wavelength  interval  ( 5) , % 


A  computer  progran  was  developed  to  calculate:  HI  which  is  the 
heat  flux  absorbed  at  the  front  surface  of  the  first  layer,  H1B  which 
is  the  heat  flux  absorbed  at  the  back  surface  of  the  first  layer 
and  H2  which  is  the  heat  flux  absorbed  at  the  front  surface  of  the 
second  layer.  Knowing  the  color  temperature  and  the  transmittance 
properties  of  pyrex,  the  equivalent  black  body  radiation  (BBR)  curve 
is  determined.  The  graphite  imaging  furnace  operates  at  a  color 
temperature  of  1800  degrees  Kelvin,  the  element  being  enclosed  in 
an  argon-filled  pyrex  bell  jar  which  limits  the  radiation  to  the 
spectral  region  between  0.29  and  3.9  microns.  By  multiplying  the 
hemispherical  black  body  radiation  curve  taken  from  the  Internation¬ 
al  Critical  Tables,  Figure  1,  by  the  transmittance  curve  of  the  pyrex 


Radiation 


f  Black  Body  Radi 


bell  jar.  Figure  2,  we  can  establish  the  black  body  radiation  dis¬ 
tribution  of  the  Incident  energy,  Figure  3.  Black  Body  Radiation 
Data,  Setting  Correction  Data,  Reflectance  and  Transmittance  Data 
and  a  copy  of  the  computer  programs  and  outputs  can  be  found  in 
Appendix  A. 

The  energy  under  the  BBR  curve  must  be  equal  to  the  incident 
energy  which  is  measured  by  a  totally  absorbing  radiometer.  This 
is  accomplished  by  multiplying  the  incremental  value  of  the  BBR 
data  by  the  incident  energy  of  the  source  and  dividing  by  the  average 
value  of  the  BBR  data,  thus  the  total  energy,  or  area  under  the  BBR 
curve,  is  equal  to  the  energy  measured  by  the  radiometer  in  confor¬ 
mance  with  the  BBR  distribution. 


The  optical  properties  of  the  two-layer  system  are  accounted 
for  by  use  of  the  equations  for  determining  HI,  H1B  and  H2: 


HI  = 


bbr(AV)N 


HI  B  = 


H2  = 


_ L 

BBR(AV)N 


A  *  3.825 
l  BBR(a)A1(a) 
A  =  0.375 


Eq.  (3) 


BBR(AV)N 


A  =  3.825 

l  BBR(a)[T1  (a )R2( a )A1  (a)  +  T1(a)R2(a)R1(a) 

R2(a)A1 (a)] 

A  =  0.375 

Eq.  (4) 

A  =  3.825 

l  BBR(a)[T1  ( a ) A2 ( a )  +  T1(a)R2(a)R1(a)A2(a)] 
a  =  0.375 

Eq.  (5) 


where 

l  refers  to  top  layer,  2  refers  to  base  layer,  B  refers  to 
back  surface 


A(/)  =  Absorptance  at  each  wavelength  interval,  % 

R(a)  =  Reflectance  at  each  wavelength  interval,  =  R( a )SC( a ) ,  % 

T(>)  =  Transmittance  at  each  wavelength  interval,  % 

>.  =  Wavelength  of  radiation,  microns 
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7©  Tronsmittonce 


Tronsmittonce  of  Pyrex  Jor 


Fia.  2.  Transmittance  of  Pyrex  Bel1  Jar 


2 

H  =  Heat  flux  perpendicular  to  surface,  cal/cm  sec 

2 

Q  =  Incident  energy  measured  with  a  radiometer,  cal/cm  sec 
BBR( a )=  Black  Body  Radiation  a4  each  wavelength  interval,  % 

88  (AV)  — J—  l  BBR(a) 

"  \  =  0.375 

N  =  Number  of  wavelengths  of  x 

Figure  4  shows  the  energy  Q  as  it  passes  through  the  two-layer 
system.  The  interface  is  separated  only  to  exhibit  the  optical  paths 
of  Q. 


where  for  all  wavelengths: 


R^a)  -  Rn(A) 

A-j  ( A )  =  Ai  i  (  a)  =  Ai g ( A ) 

T^A)  =  Tn(x)  =  T12(a) 

R2  (  a  )  =  Rg  1 

( A )  =  Agi ( A ) 

Tg( A)  ~  T2i (a) 

Where  the  first  subscript  refers  to  the  layer  and  the  second  sub¬ 
script  refers  to  the  iteration  of  reflectance ,  transmittance  or  ab- 
sorptance. 

The  values  of  HI,  and  H1B  and  H2  established  in  terms  of  the  optical 
properties  of  the  first  and  second  layer  and  the  energy  output,  color 
temperature  and  transmittance  through  pyrex  of  the  source  radiation  are 
used  in  a  Griffith  and  Horton  equation  which  considers  the  thermal  pro¬ 
perties  of  the  first  and  second  layer.  The  equation  for  calculating  tem¬ 
perature  rise  caused  by  HI,  H1B  and  H2  follows: 

U,.  (2n»l)||8/40,t 

7  n»o '  V  9  Eq- 

where 

subscript  I  refers  to  top  loyer  I,  subscript  2  refers  to  bose  loyer  2 

and  U«  Temperature  rise 

H  *  Heot  flux  perpendicuior  to  surfoce 
X  «  Totol  thichness  from  surfoce  to  point  of  temperoture 
rise  measurement 
a  *  Thichness  of  loyer  I 
D  *  Thermal  diffusely  «  h/s 
k  «  Thermal  conductivity 

S  *  Volume  specific  heot  (density  x  specific  heot) 

7  *  h|  S|  S|ht5; 
kj  S?  -  /k,  S|kt52 

A  •  (kt/D(*k,y^)*' 

t  !  Time 


Q 


Equation  6  Is  used  to  calculate  the  temperature  rise  caused 

by  HI,  H1B  and  H2.  When  calculating  the  temperature  rise  caused  by: 

HI,  x  equals  0.1032  cm  and  a  equals  0.0532  cm, 

H1B,  x  equals  0.050  cm  and  a  equals  0.0532  cm,  or 

H2,  x  equals  0.050  cm  and  a  equals  0.0  cm. 

The  total  temperature  rise  at  depth  in  the  second  layer  is 
equal  to  the  sum  of  the  individual  temperature  rises: 

U2  =  U21  +  U22  Eq. (7) 

where 

U2  =  Total  temperature  rise  at  depth  in  layer  2 

U21  =  Temperature  rise  caused  by  HI 

U22  =  Temperature  rise  caused  by  H2  plus  H1B 

H1B  is  added  to  H2  because  the  levels  of  energy  contributed  by 
H1B  were  fourH  to  be  small  in  our  particular  application. 

The  physical  and  thermal  properties  of  the  silicone  rubber  and 

the  skin  simulant  are  shown  in  Table  I. 

TABLE  1 


Physical  And  Thermal  Properties  of  Silicone  Rubber  And  Skin  Simulant 


Syri)ol 

Property 

Units 

Silicone  Rubber 
subscript  1 

Skin  Simulant 
subscript  2 

RH0 

Density 

gm/cm3 

1.0968 

1.8200 

C 

Specific 

Heat 

0.3500 

0.3571 

K 

Thermal 

Conductivity 

cal/cm°C  sec 

5.00  x  10'4 

1.31  x  10'3 

S 

Volume 

Specific 

Heat 

cal/cm3oC 

0.3839 

0.6499 

D 

Diffusivity 

2, 

cm  /sec 

13.02  x  10'4 

20.16  x  10'4 

a 

Thickness 

cm 

0.0532 

1.0 

in 


The  heat  flow  relationships  of  the  silicone  rubber  patches 


and  the  skin  simulant,  and  the  corresponding  values  for  x  and  a 
are  indicated  in  Figure  5. 

FIGURE  5 

Heat  Flow  Relationships  and  Distances 


m. 


LAYliR  1 

Silicone  Rubber 


1  AYER  2 
Skin  S'mulant 


■II1B 


XTC 


Xj  ■  distance  from  surface  where  111  is  absorbed  to  TC  *  0.1032  cm 
Bj  ■  thickness  of  layer  1  >  .0532  cm 

x j  ■  distance  from  surface  where  112  is  absorbed  to  TC  *  0.050  cm 

TC  ■  site  of  thermocouple  embedded  0.050  cm  beneath  surface  of  simulated  skin 


The  instrumentation  required  follows: 

A  Honeywell  Oscillograph  Model  1508  with  type  M40-350  gal¬ 
vanometers  is  used  to  record  experimental  data.  A  Hy-Cal  Engineering 
Model  C-1301-A-15-072  Radiometer  is  used  to  measure  the  Heat  Flux 
output  of  an  Arthur  D.  Little,  Inc.  Model  CSAF-3  Compound  Radiation- 
Imaging  System.  The  radiation  system  is  capable  of  heating  a 

2 

square  area  1/2  inch  on  a  side  with  a  flux  of  up  to  15  cal/cm  sec 
with  a  uniformity  of  +  1%.  It  consists  of  a  resistance  heated 
graphite  element  used  in  conjunction  with  a  compound  thermal 

11 


imaging  system  and  a  heat  flux-redistributor  to  provide  the 
desired  flux  uniformity  (2).  The  skin  simulant  has  a  fine  copper- 
constantan  thermocouple  embedded  .050  cm  beneath  its  surface  to 
measure  the  temperature  at  depth.  The  radiometer  and  the  skin 
simulant  thermocouple  were  calibrated  using  a  Minneapolis  Honey¬ 
well  Potentiometer  Model  2720.  Calibrations  are  made  with  a  vol¬ 
tage  standard  placed  in  series  with  each  transducer  and  the  system- 
sentitivities  are  shown  in  Table  II. 

TABLE  II 

System  Sensitivities 


Radiometer 

Skin  Simulant 

Transducer  Sensitivity 

0.4255  cal/cm2sec/MV 

23.66°C/MV 

System  Sensitivity 

0.01198  cal/cm2sec/MM 

0.3697°C/MM 

The  protocol  for  conducting  the  experiment  follows: 

1)  adjust  radiation  source  to  desired  level  of  heat  flux  output. 

2)  take  4-second  exposure  using  radiometer. 

3)  take  10-second  exposure  of  one  pigmented  patch  on  skin  simulant. 

4)  take  4-second  exposure  using  radiometer. 

5)  allow  5-minute  cooling  period  for  skin  simulant. 

6)  change  to  another  pigmented  patch. 
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7)  repeat  steps  2  through  6  until  all  pigmented  patches  are 


irradiated. 

8)  adjust  radiation  source  to  higher  level  of  heat  flux  output. 

9)  repeat  steps  2  through  8  until  completion. 

EXPERIMENTAL  RESULTS 

The  experimental  data  is  used  to  validate  two  computer  pro¬ 
grams  (Appendix  A).  Q,  the  incident  energy  measured  with  a  radi¬ 
ometer,  is  entered  into  the  H1H2  program.  Table  III  is  generated 

by  the  computer  and  gives  values  of  HI  and  H2  for  each  combination 

2 

of  pigmented  patch  and  skin  simulant  for  Q  equal  to  1.00  cal/cm  sec. 


TABLE  III 

Output  From  Computer  Program  H1H2 


W 

W+l 

W+2 

W+3 

W+4 

W+5  BLACK 

HI  = 

.34970 

.69526 

.77211 

.82475 

.84806  .87556 

H2  = 

.25776 

.16892 

.10797 

.07405 

.04530 

.03446  .02061 

HI  and  H2,  the  heat  fluxes  absorbed  at  the  first  and  second 
layer,  are  entered  into  the  DELT  program  and  the  computer  generates 
the  data  in  Table  IV,  i.e.,  time  in  seconds  vs  the  temperature  rise 
at  depth  in  °C  in  the  second  layer  caused  by  HI,  H2  plus  HI B  and  HI 
plus  H2  plus  H1B  for  the  W(white)  patch  on  the  skin  simulant  when 
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CO  — J  LTt 


2 

irradiated  with  a  heat  flux  of  1.00  cal/cm  sec.  As  observed  in 
Table  IV,  a  correction  of  -0.03  sec  is  applied  to  give  an  exact 
comparison  with  the  experimental  data  at  each  recorded  data  point 
by  compensating  for  the  shutter  opening  time  of  0.03  sec. 


TABLE  IV 

Output  From  Computer  Program  DELT 


HI  = 

.34970  H2  =  .25776 

TIME 

DELTA  T-lll 

DELTA  T-H2 

TOTAL  DELTA  T 

2.97 

5.75 

9.11 

14.86 

5.97 

13.07 

15.78 

29.65 

9.97 

22.83 

22.63 

45.46 

TIME 

DELTA  T-Hl 

DELTA  T-H2 

TOTAL  DELTA  T 

.97 

.56 

2.96 

3.52 

1.97 

2.90 

6.30 

9.20 

2.97 

5.75 

9.11 

14.86 

3.97 

8.59 

11.56 

20.16 

4.97 

11.31 

13.77 

25.07 

5.97 

13.87 

15.70 

29.65 

6.97 

16.29 

17.65 

33.94 

7.97 

18.50 

19.40 

37.98 

8.97 

20.75 

21  .06 

41.81 

9.97 

22.83 

22.63 

45.46 

Table  V  identifies  the  pigmented  patch  and  the  experimental 

and  theoretical  temperature  rises  per  unit  flux  at  3  and  6  seconds. 

2 

The  column  labelled  "Experimental  '.t/cal/cm  sec  corrected  for 
Edge  Effect"  shows  the  influence  of  edge  losses  which  is  discussed 
in  Stoll (6).  At  6  seconds'  exposure  time  a  4,  edge  effect 
must  be  made  for  the  heat  loss  at  the  site  of  measurement  because 
of  the  relatively  small  size  of  the  radiation  aperture.  This  edge 
effect  correction  increases  with  both  time  and  intensity  of  heat 


flux  exposure. 


TABLE  V 

COMPARISON  OF  EXPERIMENTAL  AND  THEORETICAL 
TEMPERATURE  RISE  DURING  IRRADIATION 


Pigmented 

Patch 

Experimental 
At/cal/cm2sec 
@  3  sec 

Theoretical 
At/cal/cm2sec 
?  3  sec 

Experimental 
At/cal/cm2sec 
@  6  sec  corrected 
for  edge  effect 

Theoretical 
At/cal/cnrsec 
@  6  sec 

HUH 

15.41 

14.86 

29.54 

29.65 

15.40 

15.08 

31.96 

32.33 

EtH  * 

15.07 

15.24 

31.84 

31.18 

19 

15.32 

15.31 

32.47 

35.15 

IxSh 

15.17 

15.16 

32.88 

35.48 

■in 

15.11 

15.16 

32.83 

35.74 

Wmsmm 

15.53 

15.12 

34.67 

35.98 

DISCUSSION  AND  CONCLUSION 


The  agreement  achieved  between  the  experimental  and  the 
theoretical  results  throughout  the  range  of  varied  optical  proper¬ 
ties  in  the  irradiated  layers  indicates  that  the  method  is  reliable 
within  experimental  limits.  The  accuracy  of  the  radiometer  used 
for  measuring  the  incident  energy  is  +  2.93%  (7)  and  the  accuracy 
of  the  galvanometer  recording  system  is  approximately  2%. 

The  fact  that  the  experimental  temperature  rise  is  always 
lower  than  the  theoretical  at  longer  exposure  times  is  ascribed 
to  a  systematic  loss  not  accounted  for  in  this  method.  Such  a 
loss  can  occur  through  scattering  or  re-radiation  from  the  surfaces. 
Much  difficulty  was  encountered  in  measuring  the  reflectance  of  the 
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white  silicone  rubber  patch.  Also,  we  suspect  that  the  thermal  and 
optical  properties  of  the  patches  change  somewhat  with  increasing 
temperature.  The  value  of  the  thermal  conductivity  constant  of 
the  silicone  rubber  was  supplied  by  the  manufacturer;  individual 
measurements  were  not  made  for  each  silicone  rubber  patch  and  it  is  ex¬ 
pected  that  each  will  vary  somewhat  from  the  value  used  in  the 
equations.  The  overall  level  of  accuracy,  about  5%, is  considered 
adequate  for  the  application  intended,  namely,  the  measurement  of 
the  thermal  conductivity  in  living  human  skin. 

In  conclusion,  the  validation  of  the  two-layer  system  for 

describing  transient  heat  flow  was  obtained  by  irradiating  a 

series  of  increasingly  dark,  black-pigmented  patches,  overlaid 

on  a  skin  simulant,  comparing  observed  and  theoretical  temperature 

rise  values.  Energy  was  supplied  by  a  graphite  imaging  furnace  and 

temperature  rise  was  measured  at  0.05  cm  beneath  the  surface  of  the 

2 

skin  simulant.  Incident  energy  levels  up  to  1.0  cal /cm  sec  were  used 
and  close  agreement  was  found  between  the  experimental  and  theore¬ 
tical  results.  Validation  of  the  equations  ensures  that  measure¬ 
ments  of  the  temperature  rise  at  the  interface  of  the  rubber  patch 
and  intact  living  skin  may  be  used  to  determine  the  thermal  con¬ 
ductivity  of  the  skin,  furthermore,  s.ince  the  primary  site  of  the 
skin  burn  injury  is  trie  epidermal-dermal  interface,  and  all  the 
pigment  is  contained  in  the  epidermis,  the  system  may  also  be  adapted 
to  evaluations  of  the  effect  of  natural  pigmentation  by  treating  the 
skin  itself  as  a  two-layer  system. 
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APPENDIX  A 


Equation  1  -  Absorptar.ee  at  each  wavelength  interval  in  layer  1. 
A1  (A)  =  1  -  [R1 (a)  +  T1(A)] 


Equation  2  -  Absorptance  at  each  wavelength  interval  in  layer  2. 

A2(a)  =  1  -  [R2( a)  +  T2( a)] 

Equation  3  -'Heat  Flux  absorbed  at  the  front  surface  of  the  first  layer. 


HI 


bbr(AV)N 


A  =  3.825 
l  bbr(a)  A1 (a) 
A  =  0.375 


Equation  4  -  Heat  Flux  absorbed  at  the  back  surface  of  the  first  layer. 
n  a  =  3.825 

H1B  =  D-DD-  *  -  tt  l  bbr(a)[T1  (a)R2(a)A1  (a)  +  T1(a)R2(a)R1(a)R2(a)A1(a)] 
BBR(AV)N  A  =  0.375 


Equation  5  -  Heat  Flux  absorbed  at  the  front  surface  of  the  second  layer. 
n  A  =  3.825 

H2  =  RRR  - N  l  BBR(a)[T1 (a)A2( a)  +  T1 ( A )  R2(  A ) R1 ( A ) A2( A ) ] 

BBR(AV)N  A  =  0.375 


where 

1  refers  to  top  layer,  2  refers  to  base  layer,  B  refers  to  back 
surface, 
and 

A(A)  =  Absorptance  at  each  wavelength  interval,  % 

R(a)  =  Reflectance  at  each  wavelength  interval,  -  R(a)SC(a),  % 

T ( A )  =  Transmittance  at  each  wavelength  interval,  % 
a  =  Wavelength  of  radiation,  microns 

SC( a ) =  Setting  correction  at  each  wavelength  interval  (5),  % 


p 

H  =  heat  flux  perpendicular  to  surface,  cal/cm  sec 

O 

Q  =  Incident  energy  measured  with  a  radiometer,  cal/cm  sec 
BBR(\)  =  Black  Body  Radiation  at  each  wavelength  interval,  % 


N  =  Number  of  wavelengths  of  a 
Equation  6  -  Temperature  Rise  at  Depth  in  Layer  2. 


U« .  2ys !  e- 1 »-  « 1 1  (2«  ♦  I  )|}*/^D«t 


-[x-a|l-y^7D,(2n*l)|]  (l-erf 


subscript  l  refers  to  top  layer  l,  subscript  2  refers  to  bose  layer  2 
ond  U  *  Temperature  nse 

H  =  Heat  flux  perpendicular  to  surface 
X  =  Total  thickness  from  surface  to  point  of  temperature 
rise  measurement 
a  *  Thickness  of  layer  I 
D  *  Thermal  diffusivity  *  k/S 
k  =  Thermol  conductivity 

S  1  Volume  specific  heat  (density  x  specific  heot) 

Y  1  k;  S;  ♦/k!  5,  k|S; 

kj>  Sj  -  /k|  5|ktS{ 

X  •  (kj/Ork./B,)" 
t  1  Time 


erf  (Z)  =  2/ / 


A-2 


Total  Temperature  Rise  at  Depth  in  the  Second  Layer 


where 

=  Total  temperature  rise  at  depth  in  layer  2 
U2l  *  Temperature  rise  caused  by  HI 
U22  =  Temperature  rise  caused  by  H2  plus  H1B 


Eq. (7) 


A-3 


The  programs  are  written  for  a  Control  Data  Corporation  6600  Computer  System  using 
Fortran  Extended  language. 


Computer  Program  to  Calculate  HI  and  HI B  plus  H2: 


100  = 

P  7 ; u  1 1 A  M  H  1 1 1 2  ( I H  P  U  T  f  m  1 T  P  U  T ,  T  A  P  !•:  1 ) 

II 

DlHOiJOU-ii  i;l.R(70)  ,:;C(70)  ,R0(  14,70)  ,SO|{(70) 

l  :■*.(  = 

riMiNKlFi;  HI  (1,10)  ,Hlb  (7,70)  ,Hi:(7,7fi) 

ll 

O 

*0 

u  1 1 : ::  m  ;  i  o  i.i  t  n  i  <  .ch  i  ( ? ) ,  t  h  i jh  i  i » ( 7  > ,  t  n  u  1 2  ( 7 ) 

1  40  = 

0 1  Iv  OOOI 01.  i((l  ,10)  ,T(7,70),A(7,70)  ,00  A  (70)  , 

150  = 

Hioior  i 

1  60  = 

1'Zt?  (1,3)  (1  d-d)  ,1  =  1,70) 

1  7  '  = 

(1,3)  (00(1)  ,1  =  1,70) 

100  = 

1  ZAP  (1,3)  <(li0(K,I),i  =  l,70),K=l  ,14) 

1  °G  = 

RGAF  (1,5)  (..MKI  )  ,1=1  ,70) 

000  = 

WJ\l  (1,3)  (1.0T(I  ),I  =  1  ,70) 

r.iC  = 

1 

Fk -i.i .‘AT  ( IX ,*'.<  =  *) 

u  LJ.  0  = 

3 

Ft) Hi -  AT  ( 1  0 ( F 5 . 4 , 1  /’ ) ) 

000  = 

p 

FORI’AT  (FI  5.0) 

,SGT(70) 


tO  (.70) 


250  = 


READ  5,G 


260  = 

SUM=0.0 

270= 

DO  7  1=1 f 70 

200  = 

SUM=SUM+BBR<I> 

290= 

7 

CONTINUE 

300  = 

B=SUM/70.0 

310= 

TCT=0.0 

320= 

D=U/Q 

330= 

DO  9  1=1,70 

340  = 

AQ (I ) =BHR (I )/D 

350  = 

TGT=TOT+AQ<I ) 

360= 

9 

CONTINUE 

370= 

C=T0T/7O.O 

380= 

PRINT  11, C 

390= 

11 

FORMAT (*  Q=*,lX,FR.6,//> 

400= 

J=0 

410= 

DO  15  K= 1 , 1 3 ,2 

420= 

J=J+1 

430= 

DO  13  1=1,70 

R(J,I)=RQ<K,I)*SC(I> 


13  CCNTI NU£ 

15  CCNTI NUZ 

J  =  0 

DC  19  K  =  2,M,2 
J  =  J+1 

DC  17  1=1,70 
T(J,I>=RcKK,I> 

17  CPNTIIJUZ 

19  CCNTI  liUE 
DC  27  J  =  1  ,7 
Di.  25  I  r  1,70 

ACJ,I)=1 .0-cn(J,I)+T<J,I>) 
I F (A ( J,  I ))::!  ,23,23 
21  A (J, I) =0.0 
23  CONTINUE 

25  CCUTIMU.*: 


630=  iiSACI)=l  ,0-(  (SSRCI  )*SC  (I  ))+SST (1 )) 

640=  29  ce»Ti  nuj: 

650=  D0  33  J=1 ,7 

660=  DO  31  1=1,70 

670=  HI (J,I )=<bBR(I )/D)*A(J,I ) 

6f.0=  Hlb(J,I)  =  ((Br»R(I)/r)*T(J,I)*SSR(I)*SC(I)*ACJ,I))+ 

690=  1  ((bBn<I)/D)*T<J,I)*i>iiR(I)*bC(I)*n(J,I)*OOR(I)*b'C(I)* 

700=  2A(J,I )) 

710=  H2(J,I)  =  (aii»n(I)/D)*T<J,I)*SbA(I))*K(BnR(I)/r>*T(J,I)* 

720=  lbOR(I)*SCa)*R(J,I)*BbACI)>  +  ((EibH(I)/D)*T(J,I)*SSR(I)* 

730=  2bC (I )*R(J, I )*OSR(I )  *bC (I )*R( J, 1 ) (I  )) 

7  40=  31  Ok/hTIIllliC 

75ii=  33  orNTIiUU 

760=  L=0 

770=  DC  41  J= 1 ,7 

700=  bllf'M  1=0.0 

790=  MJOiUi  =0.0 

soo=  i;ur-H2=o.o 

r.  l  o  =  l=l+i  A  ~1 


020  = 

D0  35  1=1,70 

C30  = 

SUMH1 =5UMH 1+H 1  (J,I ) 

G40  = 

35 

CONTINUE 

850  = 

TRUiZHl  (L)  =SUMH  1/70.0 

PG0  = 

DO  37  1=1,70 

zz 

o 

It 

SUMH 1D=CUMM 1 b+H ID ( J , I ) 

880  = 

37 

CONTINUE 

P?0  = 

TRUH  lb  (L)  =b*UHH  IB/70,0 

900  = 

DO  39  1=1,70 

910  = 

8lUlH2=S!)rH2+H2(J,I) 

920  = 

39 

CONTINUE 

930  = 

TRUEM2  (L)=  (Tl!(JH2/70*0)+TKU)ilD  (L> 

940  = 

41 

CONTINUE 

950  = 

PRINT  43 

9P>0  = 

PRINT  45,(TRUEH1(L)  ,L  =  1,7) 

970  = 

PRINT  49, (TRU£H2(L ) ,L= 1 , 7)  • 

900  = 

43 

FORMA T ( 1  OX , *W* , 7X , * V+ 1  * ,  6X , * W+2* 

990  = 

1  6X  ,'*W+4*  ,6X,  *W+5* ,  5X  ,*BLACK* ) 

1000  = 

45 

FORMAT (*  !!!=*, 7F9. 5, /> 

A-8 


1010  = 


49  FORMAT (*  H2=*,7F9.5,/> 


1020=  2ND 

Output  Generated  by  H1H2  Program: 


Q=  1.000000 

• 

w 

W+l 

W+2 

rl  1  =  .34970 

.55450 

.69526 

H2=  .25776 

•  1  6H92 

.10797 

W+3 

• 

W+  4 

W+5 

HLACK 

7721  1 

.02475 

.84206 

.07556 

07405 

.04530 

.03446 

.02061 

Computer  Program  to  Calculate  Temperature  Rise  at  Depth  in  Layer  2: 

100=  PROGRAM  D£LT( INPUT, OUTPUT) 

110=  DIf'£l2.I0N  112(100)  ,U2l  (100)  ,U22(  1 00) , TI i;£(  1  00 ) 

120=  KJAL  K 1 ,K 2, LAMBDA 

130=  i!H01  =1  .0968 

14C  =  h'HO.'l=l  .r  joo 

150=  C 1 =0.3500 

160=  02=0.3571 

170=  h'l=5.00£-4 

100=  K2= I .31 £-3 


A-9 


190  = 

XI =0.1032 

ii 

o 

o 

CJ 

ALPHA1 =0,0532 

210  = 

X2=0,05 

220  = 

ALPHA3=0,00 

230  = 

3 

F0RMATCF1 5,0) 

240  = 

5 

FORMAT (2X,*H  1  =*) 

250  = 

7 

FORMAT (  2X  ,*H2=+  ) 

ro 

3) 

o 

II 

PR  HIT  5 

270  = 

R£AC  3, HI 

220  = 

PRIM  7 

KO 

o 

II 

RrlAR  3,H2 

300  = 

SI =KHC1*C1 

310  = 

G2=RH02+C2 

32C  = 

Dl  =K 1 /S 1 

330  = 

U2=K2/S2 

340  = 

LAMM DA = 1 ,C/(K2*S0RT(D1 >-Kl*CaRT<D2>> 

350  = 

GAMf/A  =  (K2tS2+SQPvT  (K  1*S1*K2*C2) )  / 

360  = 

l  (K2*52«SQRT(K 1 *S1*K2*S2) ) 

370  = 

E=-l ,0/GAKMA 

A- 10 


380  = 


MAX  =  10 


390=  DC  11  1=1, MAX 

400=  TIMECI )  =1-0  .*^3 

<410=  SUM  1=0.0 

420=  N=-l 

430=  33  CCIITIliUE 

440=  N=N+1 

450  =  A  1  =  (2.0+H1  +  LAM13DA*  SQRT(Dl))  /GAMMA 

400  =  H  1  =2.0*SQRT  (l'2*TIME (I  )/3. 1  41  5927) 

470=  FI =(X1 )-C (ALPHA!)* Cl .0-(SQRT(D2/Dl )*<(2.C*N)+1 .0)) )) 

4F0  =  G1=-(F1**2)/(4.0*D2*TIME(I>> 

490=  El  =F1  / (2.0*$QRT (D2*TIME(I ) ) ) 

500=  P1=(G**N)*C (B1*EXP<GI))-(F1*C0MERF<E1>)> 

510=  5uru=sum+Pi 

520=  IFCN-5) 33,33,35 

530=  35  CO MIU UE 

540=  U21 (I ) =Al*SUM 1 

550=  11  CONTINUE 

5*0  =  DC  19  1  =  1, MAX 


A-n 


intmmm 


TIME(I)=I-0.03 


SUM2=0.0 


N=-l 


37  CONTINUE 


N=N+1 


A2=C2.0*H2*LAMBDA*SQRT  (D  1  )> /GAMMA 


B2=2. 0* SORT  (D2*TIMZ  (I.  )/3. 1415927) 


F2= (X2) - C  (ALPHA 3)*  ( 1  .O-CbQRT (D2/D1  >*  ( (2*  IJ)  +  1 )  )  )  ) 


G2=-<F2**2)/<4.0*D2*TIM£<I>> 


E2=F2/(2.C*i;QKT<D2*TIM£<I))> 


P2=  (0**  N>*<<fi2*£XP<G2>>-<F2*Ce»12RF<E2>>> 


SUM2=SUM2+P2 


IF(N-5)37, 37,39 


39  CONTINUE 


U22(I)=A2*GUM2 


19  CCNTINUE 


PC  21  1  =  1, MAX 


U2 ( I )  =  U2 1 ( I ) +U22 ( I ) 


5 


21  CONTINUE 


PRINT  23,H1,H2 


23  FKRMAT(2X,*Hl=*,F8.5,2X,*H2=*,F8.5,/> 


PRINT  25 


25  F0MAT<3X,*TIME*(HX.*D£LTA  T-Hl*  ,8X,*DELTA  T-H2*,8Xt 


1+T2TAL  DELTA  T*,/> 


pr(  i  (ji  27  ,  TIN  EC  3)  ,1121  (3),U22(3)  ,U2(3) 


PRINT  27,Tir£(6),U2l(6) ,U22(6) ,U2(6) 

PRINT  27 ,TINE( 1 0) ,U2 1 ( 1 0) , U22( l 0) ,U2 ( l 0) 

27  format ( 2x ,F5 . 2, 1  lx , fs .2 , 1  ax ,F6. 2, 1 2X ,F6 . 2 > 


PRINT  29 


29  FCRt1AT(/> 


PRINT  25 


D2  31  1  =  1 /‘AX 


PR 


I  NT  27,TIMZ(I) ,U2l (I ) ,U22(I ) ,U2(I > 


31  CONTINUE 


PRINT  25 


END 


FUNCTI3N  CCMERF  <U) 


APPRKXinATIGIl  FORMULA  “  CCP.PLE'' ENTARY  t RR&P-  FUNCTIf 


•■o  ut  r- 


940=  Ct)M£RF=  1  ./((((  C  (  .430633£-4*U+2.765672E-4)*ll+l  .  52  A  43 E-4) 

950=  l*U+92.705272£-4)*U+422.820l23£-4>*ll+705.2307n4£-4>*IRl  .)* 

960=  2*16. 

970=  RETURN 

920=  END 


Output  Generated  By  DELT  Program: 

W(White)  Silicone  Rubber  Patch  on  Skin  Simulant 


H 1 =  .34970  H2  =  .25776 

TIME  DELTA  T-Hl  DELTA  T-H2 


.97 

.97 


5.75  9.11 
13.27  15.72 
22.23  22.63 


TIME 

DELTA  T 

.97 

.56 

1  .97 

2.90 

2.97 

5.75 

3.97 

8.59 

4.97 

1  1  .31 

5.97 

13.07 

6.97 

16.29 

7.97 

12.50 

11.97 

20.75 

9.97 

oo  r»  T 

•  o  O 

HI  DELTA  T-H2 

2.96 
6.30 
9.11 
11.56 
13.77 
15.72 
17.65 
1  9 . 40 
21  .06 
22.6  3 


TOTAL  DELTA  T 

1  4.D6 
29.65 
45.46 


TOTAL  DELTA  T 

3.52 
9.20 
14.06 
20.16 
25.07 
29.65 
33.94 
37.92 
41  .31 
45.46 
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W  +  1  Silicone  Rubber  Patch  on  Skin  Simulant 


HI:  .55450  M?.=  .  1 


II  HE 

D£LTA  T-Hl 

DELTA  T-Hid 

TDTAL  D £L T A  T 

2.97 

9.11 

5.97 

15. on 

5.57 

121  .99 

10.34 

312.33 

9.57 

36.120 

1  4.R3 

51  .03 

Tine: 

D£LTA  T-H  1 

DELTA  T-Hli 

TOTAL  DELTA  T 

.57 

.co 

1.94 

2.  n  3 

1  .97 

4.59 

4.13 

r  .72 

12.97 

9.11 

5.97 

1  5.09 

3.97 

1  3 . 63 

7.50 

2 1  . 20 

4.57 

17.93 

9.02 

26.95 

5.97 

:;i  .99 

10.34 

32.33 

6.97 

125.012 

11.57 

37.39 

7.97 

119.45 

1 12 .  72 

42.17 

f.97 

32.91 

i3.no 

46.71 

5.  "'7 

36.no 

1  4.23 

51  .03 

W  +  2  Silicone  Rubber  Patch  on  Skin  Simulant 


1=  .65526 

illl:  .10797 

TI!  £ 

DELTA  T-SU 

D£L  FA  T-H 2 

TOTAL  D£LTA  T 

2.97 

11.43 

3.02 

15.24 

3  #3  7 

127.57 

6.61 

34,10 

3.97 

45.39 

9 . 40 

54.07 

T I  i 1  „ 

D£LTA  T-H 1 

DELTA  T-H2 

TOTAL  DELTA  T 

.97 

1.11 

1.24 

2.35 

1.97 

5.76 

2.6,4 

0.40 

2.97 

1  1  .43 

3.02 

1  5.24 

3.97 

1  7 . 00 

4.04 

21.93 

4.97 

O  ^  • 

c-  4.-  • 

5.77 

t  >  '*•  ri  c 

5.97 

27.57 

6.61 

3  -t .  1  Q 

6.97 

32.32 

7.39 

39.77 

7.47 

36.93 

8.13 

*<5.06 

0.97 

41.26 

2.02 

c  r.  r  .  ^ 
u  •  wi  i 

9. 97 

M  >  •  vj  J 

9.  40 

34.07 
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W  +  3  Silicone  Rubber  Patch  on  Skin  Simulant 


Hl=  .7721 1  H2:  .07405 


Tint 

DELTA  T-Hl 

DELTA  T-H2 

TflTAL  DELTA  T 

2.97 

12.69 

2.62 

15.31 

5.97 

30.62 

4.53 

35.15 

9.97 

50.41 

6.50 

56.91 

TIME 

DELTA  T-Hl 

DELTA  T-H2 

TOTAL  DELTA  T 

.97 

1  .23 

•  G5 

2.00 

1  .97 

6.39 

1.01 

0.21 

2.97 

12.69 

2.62 

15.31 

3.97 

12.97 

3.32 

oil  QQ 

IjL.  .  . ,  J 

4.97 

24.96 

3.95 

20.92 

5.97 

30.62 

4.53 

35.15 

6.97 

35.96 

\  5.07 

41 .03 

7.97 

41  .01 

5.57 

46.59 

0>.  97 

45.22 

6.05 

51.07 

9.97 

50.41 

6.50 

56.91 

W  +  4  Silicone  Rubber  Patch  on  Skin  Simulant 


HI:  .22475  H2=  .04530 


TIKE 

DELTA  T-Hl 

DELTA  T-H2 

TOTAL  DELTA 

£.97 

13.56 

1.60 

1  5.16 

5.97 

32.71 

2.77 

35.40 

9.97 

53.04 

3 . 90 

57.02 

TILE 

DELTA  T-Hl 

DELTA  T-H2 

TCTAL  DELTA 

.97 

1  .  7  2 

.52 

1  .04 

1  .97 

6.03 

r.u 

7.9  4 

2.97 

13.56 

1  .60 

15.16 

3. 97 

20.27 

o  n 7, 

22.30 

4.97 

26.67 

iL  • 

.29.00 

5.97 

32.71 

2.77  . 

35.40 

6.97 

30.41 

3.10 

41.51 

7.97 

43  •  y  l 

3.41 

47.22 

0.97 

40.9  4 

3.70 

52.65 

9.97 

53.04 

3 . 93 

57.02 
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W  +  5  Silicone  Rubber  Patch  on  Skin  Simulant 


Hl=  .04806  H2=  .03446 


Tint 

DELTA  T-H 1 

DELTA  T-H2 

TOTAL  DELTA  T 

2.37 

13.94 

1  .22 

15.16 

6.97 

33.63 

2.11 

35.74 

3.97 

55.37 

3.02 

50.39 

Tir-E 

DELTA  T-H  1 

DELTA  T-H 2 

TOTAL  DELTA 

.97 

1  .35 

.40 

1 .75 

1  .97 

7.02 

.34 

7.27 

2.97 

13.94 

1  .22 

15.16 

0,9  i 

20.3  4 

1  .55 

22.39 

*4.37 

27.42 

1  .34 

29.26 

5.97 

33.63 

2.11 

35.74 

6.97 

39.49 

2.36 

41  .35 

7.97 

45.05 

2.59 

47.64 

2.97 

30.33 

L«C2 

53.14 

9.97 

55.3  7 

3*02 

52.3  9 

Black  Silicone  Rubber  Patch  on  Skin  Simulant 


Hit  .07556  H2  =  .02061 


TIDE 

DELTA  T-H  1 

DELTA  T-H2 

TOTAL  DELTA 

2.97 

14.39 

.73 

1  5.12 

5.97 

34.72 

1.26 

35.93 

9.97 

57.16 

1  .31 

58.97 

TIME 

DELTA  T-Hl 

DELTA  T-H2 

TOTAL  DELTA 

.97 

1  .40 

.24 

1  .63 

1  .97 

7.25 

.50 

7.76 

2.97 

1  4.39 

.73 

15.12 

3.97 

21  .52 

.92 

22.44 

4.57 

28.31 

1.10 

29.41 

5.97 

3  4.7  2 

1  .26 

35.92 

6.97 

40.78 

1  .41 

42.19 

7.97 

46.51 

1  .55 

48.06 

3.97 

51  .96 

1  .68 

53.64 

5.97 

57.16 

1  .81 

58.97 

A- 17 


Data  Is  presented  In  groups  of  70  data  points  corresponding 
to  70  wavelenghts  from  0.375  to  3.825  microns  as  Illustrated  below: 

10-0.375,0.424,0.475,0.525,0.575,0.625,0.675,0.725,0.775,0.825 

20-0.875,0.925,0.975,1 .025,1 .075,1.125,1 .175,1 .225,1 .275,1 .325, 

30=1.375,1.425,1.475,1.525,1.575,1.625,1.675,1.725,1.775,1.825, 

40-1.875,1.925,1.975,2.025,2.075,2.125,2.175,2.225,2.275,2.235, 

50-2.375,2.425,2.475,2.525,2.575,2.625,2.675,2.725,2.775,2.825, 

60-2.875,2.925,2.975,3.025,3.075,3.125,3.175,3.225,3.275,3.325, 

70-3.375,3.425,3.475,3.525,3.575,3.625,3.675,3.725,3.775,3.825, 

The  Data  File  Key  used  In  the  H1H2  program  follows.  Each 
data  point  has  a  corresponding  wavelength  associated  with  It  as 
described  previously. 

Data  File  Key 

Lines  100  to  160  -  Black  Body  Radiation  Data  -  BBR(x) 

Lines  180  to  240  -  Setting  Correction  Data  -  SC(x) 

Lines  260  to  320  -  Reflectance  of  White  Patch  -  R(x) 

Lines  330  to  390  -  Transmittance  of  White  Patch  -  T(x) 

Lines  400  to  460  -  Reflectance  of  W  +  1  Patch  -  R(x) 

Lines  470  to  530  -  Transmittance  of  W  +  1  Patch  -  T(x) 

Lines  540  to  600  -  Reflectance  of  W  +  2  Patch  -  R(x) 

Lines  610  to  670  -  Transmittance  of  W  2  Patch  -  T(x)  • 

Lines  680  to  740  -  Reflectance  of  W  +  3  Patch  -  R(x) 

Lines  750  to  810  -  Transmittance  of  W  +  3  Patch  -  T(x) 

Lines  820  to  880  -  Reflectance  of  W  +  4  Patch  -  R(x) 
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Lines  890  to  950  -  Transmittance  of  W  +  4  Patch  -  T<x) 

Lines  960  to  1020  -  Reflectance  of  W  +  5  Patch  -  R(x) 

Lines  1030  to  1090  -  Transmittance  of  W  +  5  Patch  -  T(x) 

Lines  1100  to  1160  -  Reflectance  of  Black  Patch  -  R(x) 

Lines  1170  to  1230  -  Transmittance  of  Black  Patch  -  T(x) 

Lines  1250  to  1310  -  Reflectance  of  Skin  Simulant  -  SSR(x) 

Lines  1330  to  1390  -  Transmittance  of  Skin  Simulant  -  SST(x) 


TAPE  1 


100=.0002,.0007,.0030,.0090,.0214,.0423,.0817, . 1 1 64, . 1 679, .2297, 
1 10=. 2944,. 3620,. 4302,. 4834,. 5243,. 5645,. 6089,. 651  5,. 6942,. 7372, 
1 20= . 7750,. 8209,. 8691 ,.8878,. 8833,. 8758,. 8682,. 8607,. C541 ,.8446, 
130=. 8380,  .8368, .8227, .8062, .7837, .751 1 , .721 6, .6982, .6883, .6664, 
140=. 6300,. 5981 ,.5718,. 5451 ,.5267,. 4946,. 4303,. 0624,. 01 89,. 0334, 
1 50=. 0523,. 0761,. 088 4, .1079,. 11 94,. 1298,.  1  423,.  1467,. 1397,. 1  186, 
160=. 0729,. 0320,, 0072,. 0031 ,.0000,. 0000,. 0000,. 0000,. 0029,. 0025, 
1 70=*ECR 


180=. 97 63,. 9688 ,.96 50 ,.96 50, .9625, .9575, ,953 7, .9 51 3, .9475 
190=. 9 400,.  9400,. 9372,. 93 18, .9277, .9252, .9228, .9202, .9185 
200= .91 65, ,91 55, .91 38, .9112,. 9093, .9078, ,9073 , • 9078 , . 9073 
21 C=. 9020,  .8960, .8920,. 8900, .8863 , .8807, .8762, .8727 , .8690 
220=. 8645, .8675, .8682, .8667, .8690, .8750 , .88 1 0, .8800 , .8800 
230= .8800,  .8800, .8800, .3800, .8800, .2800, .8800, .3800, .8800 
240  =  . 8800,.  8800,. 8800,. 8800,. 8800,. 8800,. 8800, .8800,. 8800 


,.9425 
,.9175 
,  .9058 
,  .8650 
,.8800 
, .8800 
, .8300 


t 


» 


t 


t 


25O=*£0R 

260  =  . 9999,.  9999,. 8554,. 8293,. 7698,. 7578,. 7090,. 681  9, .64  54, .631 2, 
270=. 6210, .5887, .5701  ,. 5265, . 5078 ,. 4742, .381  1  , .423  5, .4203, .3861 , 
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280= .2838,  .2316, .2966, .2443, .2921 , .2791 1  765,  .0634, . 1  1 70, . 1 683, 
290=. 21  06,.  1 664,. 1 438,. 1 401 ,.1516,. 1398,. 0626,. 0000,. 0000,. 0000, 
300= .0000,  .0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000, 
310=. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000, 
320=. 0000,  .0000, .0000,  .0000, .0000, .0000, .0000,  .0000, .0000, .0000, 
330=. 2510,. 3360,. 3570,. 3680,. 37 40,. 3890,. 41 20,. 4270,. 4430,.  4550, 
340=. 4640, .4770, .491  0,  . 5090 ,. 5260, . 5280, . 4690, .  5490, . 5600, .5670, 
350=. 51 00,. 51 60,. 5610,.  5490,. 5960,. 6020,. 4720,.  41  20,. 5280,. 6000, 
360=.  6460,  .6610,  .6520,  .6620,  .6830,  .6600, .  568,0,  .  4290,  .  2400 ,  .3  1  90, 
370=. 3320,. 3890,. 3930,. 0510,. 0600,. 03 50,. 0250,.  01  90,. 0000,. 0000, 
380=. 0000, .0000, . 0000,. COOO,. 0000,. 0000, .0000, .0000, .0000, .0000, 
390=. 00 00, .0000, .00 00, .0080, .0160, ,08 50, .14 00, .1400, ,1500, ,8 130, 
400=. 23 10,, 2330,. 2260,. 21  10, ,2010, . 1940, . I  920, .1900,. 1890, ,1  920, 
4 10=. 2080,. 2460,. 2890,  .2950, .271 0,.2380, .2030,  .2080, .21 20,. 2280, 
420= .2380, .2610, .2950, .2990, . 3 1 30, . 3090, .2420, .2020, .2380,  .2580, 
430= .2640, .2560 , .2400, . 2330, . 2280 ,.20 90, .1690,  .11  90, .0710, .0810, 
440= .0810,  .0860, .0830, .0900, .0800, .0590, .0590,  .0590, .0590, .0590, 
450=. 0000,. OOOC,. 0000,  .0000,. 0000,. 0000,. OOOC, .0000, .0000,  .0000, 
460=. 0000, .0000, .0000,  .0000, .0000, .0000, .0000,  .0000, .0000, .0000, 
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470=. 0280, .0370,. 0400,. 0430,. 0470,. 0500,. 0560,. 0610 
430= . 0930, • 1310, • 1 780, .2010, .  1  900, • 1 720, . 1 520, . 1*90 
490=. 2280, .2630,. 31 20,. 3270,. 3610,. 3700,. 3000,. 2690 
500=. 3880,. 3960,. 3930,. 4000,. 4090,. 4000,. 3530,. 2770 
510=. 21 90, .2580,. 2640,. 3 130,. 3030,. 2 100,. 21 00,. 0230 
520=. 0000, .0000,. 0000, .0000, .0000, .0000, .0000, .0000 
530=. 0000, .0000,. 0000, .0000, .0000, .0000, .0000, .0000 
540=. 1 530,, 1 51 0,. 1470,.! 370,. 128Q,. 1220,. 1200,. 11 80 
550=.  1360, .1700,, 21 20,. 21 70,.  1920,. 1630,. 1400,. 1390 
560=, 1760, .2030, .2330, .2420, ,2570, .2520, .2080 , .1770 
570=. 2080, .1990,. I860,. 1790,. 1720,. 1590,. 1370,. 1010 
530=. 07 50, .0790, .07 70, .08 10, .0760, .0530, .0580, .0580 
590=. 0000, .0000, .0000, ,0000,. 0000, .0000, .0000,. 0000 
600=. 0000,. 0000,. 0000, .0000,. 0000,. 0000,. 0000,. 0000 
610=. 0040, .0060, .0030,. 0080,. 01 00,. 01 00, .01 20,. 01 40 
620=. 03 10, .0490,. 0800,. 0960,. 0860,. 0710,. 0620,. 0680 
630=. 1  1 50, .1  480, .1860, .2040, .2320, .2390,  .  1 9R0, . 1 780 
640=. 2550, .2600, .2580,. 2620,. 2680,. 2620,. 2370,. 1820 
650=. 1 480, .1740,. 1810,. 21CC,. 2100,.! 440,. 1440,. 0090 


,.0690, .0790, 
,.1830, .21 10, 
, .3330, .3670, 
,.1520, .2030, 
, .02SO, .0280, 
, .0000, .0000, 
, .0000, .0000, 
,.1170, .1210, 
,.l 420, .1 600, 
,.1980, .2060, 
, . 06P0 ,  .0761), 
,.0580, .0530, 
, .0GC0, .0000, 
, .0000, .0000, 
,  .0180, .0210, 
,.0770, .0960, 
,.2190, .2400, 
,.1010, .1410, 
, .0090, .0090, 
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660=. 0000,. 0000,. 0000,. 0000, ,0000,. 0000,. 0000,. 0000,. 0000, 
670=. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000, 
680=.1210,.l 1  GO,. 1  120,. 1030,. 09G0,. 0930,. 0900,. OBBO,. 0880, 
690=. 1020, .13 50, .1710, .17 60, .1510, .1230, ,1070, .10 40, .1090, 
700=, 1420, .1600, .I960, .2070, .2170, .2110, ,1760, .1530, .1660, 
7 1 0=. 1720,. 1640,. 1540,. 1480,. 1420,. 131 0, . 1 1 20, .0B80, . 0630, 
720=. 0700, .0730, .0710, .0720, .0740, .0580, .0580, .0580, .0580, 
730=. 0000, .0000,. 0000,. 0000, .0000,. 0000,. 0000, .0000,. 0000, 
740=. 0000, .0000, .0000, .0000, .0000, .0000,. 0000,. 0000, .0000, 
700=. 0000, .0000, .0000, .0000, .00 00, ,00 00, .0020, .00 40, .0060, 
760=. 01 20, .0230,. 0430,. 0530, .0400, .0390, .031 0, .0340, .0400, 
770=, 0680, ,0920, .1220, .1400, ,1610, .1600, .1390, .12 50, ,1530, 
7B0=, 1  7 90, .10 10, .1800, .1820, .1880, .1820, .1630, .1280, .0700, 
790=. 1030,  .1230,. 1280,. 1550,. 1490,. 1010,. 1010, .0000,. 0000, 
800=. 0000,  .0000, .0000, .0000, .0000, .0000, .0000, .0000,. 0000, 
CIO:. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000, .0000, 
820=. 1030,  .1010,. 0980,. 0900,. 0860,. 0810,. 0790,. 0770,. 0770, 
830= . 0890,  .1 130,. 1490,. 1 540,. 1380,.!  100,. 0970,. 09 40,. 0980, 
f 40= . I  290,. 1  520,.  1790,.! 900,. 1  980,. 1930,. 1680,. 1480,. 1  550, 
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.0000, 

.0000, 

.0910, 

.1230, 

.1720, 

.0710, 

.0530, 

.0000, 

.0000, 

.0090, 

.0520, 

.1690, 

.0970, 

.0000, 

.0000, 

.0000, 

.0790, 

.1110, 

.1570, 


850=. 1550,. 1490,. 1410,. 1340,. 1280,.  1200,. 1060,. 0880,. 0670,. 0700, 
860=. 0710,. 0720, .0700, .0720, .071  0,  .0570, .0570, .0570, .0570, .0570, 
870  =  . 0000, .0000, .0000, .0000, .0000,  .0000, .0000, .0000, .0000, .0000, 
080= .0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000, 
890=. 0000,. 0000, .0000,. 0000,. 0000,  .0000, .0000,. 0000, ,0000,. 0020, 
900=. 0040,. 01 00, .0200, .0270, .0210, .0180, .0130, .01 40, .0180, .0260, 
910=, 03 70, . 053 0,. 07 40, . 0880, .1  020,  .1070, .08 90, ,08 10, ,09 90, .1080, 
920= .11 20, .11 20, .1110, .11 20, .11  40, .1110, .1 000, .0770, .0410, .0590, 
930=. 0620,. 0750, .0790,. 0980,. 0930,.  0620,. 0620,. 0000,. 0000,. 0000, 
940=. 0000, .0000,  .0000,. 0000,. 0000,  .0000, .0000,. 0000,. 0000, .0000, 
950=. 0000,.  0000,. 0000,. 0000,.  0000,.  0000,. 0000,. 0000,. 0000,. 0000, 
960=. 0940,. 0910,. 0880,. 0800,. 0760,. 0720, .0700, .0690, .0690,  .0720, 
970  =  . 081 0,. 1060,.  1390,.! 450,. 1230,.  1010,. 0880,. 08 50,. 0900,.  1010, 
98 0=, 1190, .1  430, .1680, .1800 ,.1860, .1820, .1580, .1430, .1480, .1470, 
990=. 1 460,. 1 390,. 1 300,. 1 230,. 1  1 90,  .  t  1  00 ,. 0990, .0840, .0630 ,  .0700, 
1000=. 0700,. 0710, .0700,. 0710,. 07 10,. 0560, .0560,. 0560,. 0560, .0560, 
1010  =  . 0000,. 0000, .0000,. 0000,. 0000,.  0000,. 0000,. 0000,. 0000,. 0000, 
1020=. 0000, .0000,  .0000, . 0000,. 0000,. 000C, .0000,. 0000, .0000,  .0000, 
1 030=. 0000, .0000, .0000, .OOOC, .0000,  .0000, .0000, .0000, .0000, .0000, 
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1 


10 40  =  .00 20, .0060, .0110, .0170, .01  20, 
1050=. 0220, .0360, .0520, .0630,. 0770, 
1  060=.  0070,  .0880, .  OB  70,.  OP.  80,.  08  90, 
1070=. 0500, .0600, .0620, .0780, .0730, 
1080=. 0000, .0000,. 0000,. 0000,. 0000, 
1 090=. 0000, .0000, .0000,. 0000, .0000, 
1 1 00= , OP  1  0 ,. 0790, ,0750 , .0690, . 0640 , 
1  1 10 =.07 20, .0^40, ,1230, .1290, .1100, 
1 1 20= . 1  1  00, . 1 330, . 1 550, . 1 660, . 1 730, 
1 1 30= . I  370,  1 290, . 1 21 0, . 1  1 50, . 1  1 00, 
I  1 40  = , 0620 , . 063C , ,0690, .0710, ,0700, 
1  1 50= .0003, .0000, .0000, .0000, .0000, 
1  1  60= . 0000, .0000, .0000,. 0000, .0000, 
1  170=. 0000, .0003, .0000,. 0000, .0000, 
1  ISOs. 0020, .0030,. 0060,. 0080,. 0060, 
11 90=. 01  00, .01  PC, .0300, .0390, .0470, 
1200=. 0530, .0530,, 0523, .0520,. 0530, 
1213=. 0300, .0370, .0393, .0480, .0460, 
1223=. 0003, .0000, .0000, .0000, .0000, 


.0100,.OORO, .0090, .0100, .0130, 
.0800,. 0670, .0610, .0730, .0810, 
.0870,. 0780,.  0590,. 03 20, .0480, 
.0500,. 0500, .0000, .0000, .0000, 
.0000,, 0000, .0000, .0000, .0000, 
.0000, .0000,  .0000, .0000, .0000, 
.0620,. 0600,. 0600, .0600, .0620, 
.08 90,. 0780, .0770, .0800, .0930, 
.  1710,. 1500,.  1380,. 1400, . 1350, 
.1030,. 0950,. 0800,. 0610, .0670, 
.0520,. 0580,. 0520, .0520, .0520, 
.0000, .0000, .0000, .0000, .0000, 
• 0000,.OOOOr.OGOO,.COOO,.OGOO, 
.0000, .0000, .0000, .0000, .0000, 
.0030,. 0030, .0030, .0040, .0070, 
.0490,. 0400,. 0380,. 0470, .0510, 
.0520, .0460,. 03 50, .0190, .0280, 
.0290, .0290,.  0000, .0000, .0000, 
.0000, .0000, .0000, .0000, .0000, 


A-25 


1230=. 0000, .0000, .0000, .0000,  .0000, .0000, .0000, .0000, .0000, 
1 240=*£CR 

1250=.  1710,. 2700,  .2650,  .3000, .3680, .4250, .4390,  .4390, .4350, 
1260=. 4200,. 4000,. 3910,. 3750,. 3710,. 3520,. 3000,. 28 50,. 2830, 
1270=. 2090,. 1390, .0390, .081 0,  .0950, .  1  090, . 1  1 30, .0950, .0960, 
1220=. 0870,. 0550,. 0320,. 0280,. 0280,. 0260,. 0230,. 0220,. 0210, 
1290=. 0210,. 01 90, .0180, .0180, .0180, .0190, .0190, .01 80, .0130, 
1300=. 0000,. 0000,. 0000,. 0000,.  0000,. 0000,. 0000,. 0000,. 0000, 
1310=. 0000, .0000, .0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000, 
l 320=* £8 H 

1330=. 0000, .0540, .1 160,. 1800,. 2460,. 2870,. 3070,. 3240,. 3390, 
1340=. 3670,. 3800,. 3920,. 4000,. 4200,. 4310,. 4290,. 4390,. 4530, 
1350=. 4480,. 3900,. 3300,. 3590,. 3980,. 4520,. 4500,. 4360,. 4560, 
1360  =  .  4560,. 2850,  .2920,. 2180,. 21  90,.  2240,. 21  40, . 1 930, . 1 690, 
1370=. 2320, . 2000,. I  240,. 1  1 40,.  I  550,.  1990,. 1750,.!  140,. 1  140, 
1 380=. 0000,. 0000, .0000,. 0000,. 0000,. 0000,. 0000,. 0000,. 0000, 
1  390= ,0000, .0000, .0000,  ,0000, .0000,. 0000, .0000,. 0000, .0000, 
1  4OO=*E0R 

1  41  0=*E?F 


.0000, 

.4300, 

.2640, 

.0990, 

.0220, 

.0180, 

.0000, 

.0000, 

.3550, 

.4640, 

.4830, 

.2040, 

.1140, 

.0000, 

.0000, 
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